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I.  INTRODUCTION 


Images  obtained  using  coherent  Illumination  contain  a  multiplicative  noise 
component  called  speckle.  The  noise  component  Is  a  result  of  Interference 
between  wavelets  scattered  from  a  surface  whose  roughness  Is  on  the  order  of 
the  wavelength  of  the  Illumination.  The  statistical  properties  of  this  noise 
have  been  developed  [1]  and  the  phenomenon  and  Its  suppression  have  o  en 
extensively  studied  [2]. 

The  presence  of  speckle  noise  in  an  imaging  system  reduces  the  resolution 
capability  of  the  system  by  as  much  as  a  factor  of  seven  [3]  and  Is  particu¬ 
larly  severe  in  imagery  of  low  contrasts  [4] .  Thus  the  suppression  of  speckle 
noise  Is  an  Important  component  of  a  coherent  imaging  system  design.  The 
design  of  suppression  techniques  can  best  be  optimized  by  basing  them  on  a 
detection  theory.  An  excellent  account  of  the  detection  theory  for  small  high 
contrast  objects  In  speckle  noise  has  been  given  by  Dainty  [5]. 

For  Images  In  speckle,  the  resolution  can  be  expected  to  improve  if  inde¬ 
pendent  Intensity  patterns,  or  looks,  are  superimposed.  This  Is  due  to  the 
fact  that  the  standard  deviation  of  the  mean  of  N-lndependent  measurement  sets 
is  given  by  a  =  ai/(N)72  where  oj  Is  the  standard  deviation  of  a  single  mea¬ 
surement  set.  The  noise  can  also  be  reduced  by  averaging  the  fluctuations 
over  some  area  A  of  uniform  Intensity  In  the  image.  If  s  Is  linear  dimension 
of  a  speckle  In  the  Image  plane,  then  the  standard  deviation  should  be  reduced 
by  a  factor  proportional  to  s/A'^. 

The  procedures  used  to  produce  the  final  speckle  smoothed  image  can  be 
divided  into  four  categories:  1)  reduction  in  the  temporal  coherence  of  the 
illumination,  2)  reduction  in  the  spatial  coherence  of  the  Illumination,  3) 
time  averaging  while  moving  the  receiver  or  transmitter  aperture,  and  4) 
integrating  adjacent  pixels.  These  techniques  trade  off  system  bandwidth, 
recording  time,  or  system  resolution  for  signal-to-nolse  Improvement. 

An  appealing  option  to  the  above  "N-look"  methods  of  speckle  smoothing  is 
to  perform  some  type  of  digital  processing  technique  in  addition  to  or  instead 
of  N-look  smoothing.  Conventional  imaging  processing  techniques  have  been 
developed  for  additive  noise;  however,  if  the  log  of  the  intensity  Is  first 
obtained,  then  speckle  noise  becomes  an  additive  noise  [6]  and  conventional 
techniques  can  be  applied  [7]. 

In  this  paper  a  digital  simulation  of  an  Image  with  speckle  noise  is 
developed  and  verified  by  comparison  to  experimentally  obtained  speckle 
Imagery.  The  digital  simulation  is  then  used  to  obtain  detection  probability 
statistics  and  to  evaluate  various  linear  and  nonlinear  digital  filtering 
techniques  for  speckle  reduction. 
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TI.  SPECKLE  THEORY 


Radiation  that  has  at  least  a  partial  degree  of  spatial  and  temporal 
coherence  will  exhibit  speckle  noise  due  to  volume  scattering  while  propa¬ 
gating  through  a  medium  such  as  the  atmosphere  or  due  to  scattering  from  a 
rough  surface-  The  physical  origin  of  this  noise  is  in  the  interference  of 
many  coherent  wavelets  that,  because  of  the  scattering,  have  traveled  over 
different  optical  path  lengths- 

In  free  space  propagation  it  is  easy  to  see  that  many  wavelets  will  con¬ 
tribute  to  the  intensity  at  an  observation  point-  In  a  perfect  imaging  sys¬ 
tem,  diffraction  causes  the  intensity  at  an  image  point  to  result  from  many 
wavelets  - 

To  obtain  speckle  we  require  only  that  a  large  number  of  scattering 
events  occur  within  the  point  spread  function  of  the  imaging  system  [1]- 

The  electric  field  at  any  point  (x,y)  in  our  observation  plane  is,  for  a 
monochromatic  wave 


E  =  a(x,y)e^'^^ 


where  a(x,y)  Is  the  complex  amplitude  of  the  field-  For  both  free  space  pro¬ 
pagation  and  imaging,  the  amplitude  is  given  by  the  sum  of  amplitudes  of  wave¬ 
lets 


N 

a(x,y)  =  I  I  Ok  I  exp  (l(}.ic)  •  (1) 

k=l 


To  obtain  the  statistics  of  the  electric  field,  we  must  solve  the  classi¬ 
cal  random  walk  problem  [1,2]-  In  general,  the  statistics  depend  upon  the 
coherence  of  the  radiation  and  the  statistics  of  the  scattering  surface  or 
medium-  For  simplicity,  we  will  assume  that  the  radiation  Is  coherent  and 
linearly  polarized-  We  will  only  consider  speckle  caused  by  scattering  from  a 
rough  surface  and  will  make  those  assumptions  necessary  to  allow  us  to  ignore 
the  details  of  the  scattering  process  [3]-  The  resulting  statistics  apply  to 
speckle  fields  that  were  generated  In  this  study- 

The  assumptions  to  be  made  are: 

1-  The  scattering  does  not  affect  the  polarization-  To  produce  speckle 
in  the  laboratory,  we  rely  upon  multiple  scattering  in  opal  glass,  which  depo¬ 
larizes  the  light-  A  polarizer  is  used  in  the  recording  of  the  speckle  field 
to  allow  us  to  retain  this  assumption- 

2-  Ok  4’k  statistically  Independent  of  each  other  and  of  all 
other  a*s  and  (J>'s-  This  assumption  allows  us  to  neglect  the  actual  scattering 
event  - 


3-  i)k  Is  uniformly  distributed  In  the  Interval  -it  to  ir-  This  assumption 
requires  that  the  scattering  surface  be  rough  compared  to  a  wavelength- 


2 


1 


4.  We  will  assume  that  the  observation  plane  Is  In  the  far  field.  The 
statistics  in  the  image  plane  are  the  same  as  the  far  field,  provided  there 
are  a  large  number  of  Independent  scattering  centers  within  the  point  spread 
function  of  the  imaging  system  [1,2].  With  these  assumptions  and  for  a  value 
of  N  large  enough  to  allow  application  of  the  central  limit  theorem,  we  can 
write  speckle  as  a  complex  Gaussian  field 

a(x,y)  =  aR(x,y)  +  i  ai(x,y)  (2) 

where  aR,ai  are  zero  mean  Gaussian  random  variables  (for  each  x,y)  with  vari¬ 
ance 

The  speckle  intensity  is  then,  simply 

2  2 

s  =  s(x,y)  =  a  (x,y)  +  a  (x,y) 

R  I 

and  phase 

,  ,  _i  ai(x,y) 

9  =  0(x,y)  =  tan  ^  ^ 

aR(x,y) 

The  random  variables  9  and  s  are  found  to  be  Independent  with  densities 

p(s)  =  — ^  exp  (-s/2a^)  ,  s  >  0  ,  (5a) 

2a 


(3) 


(4a) 


(4b) 


p(9)  =  V2Tr,  -n  <  9  <  IT  , 

(5b) 

p(s, 9)  =  p(s)  p(9)  . 

(5c) 

The  intensity  at  a  point  in  the  observation  plane 
tlal  statistics  and  the  phase  obeys  uniform  statistics, 
and  phase  are  statistically  independent.  We  will  call 
speckle . 

obeys  negative  exponen 
.  Also  the  Intensity 
this  fully-developed 

The  Intensity  is  the  parameter  measured  experimentally  and  thus  of  pri¬ 
mary  Interest.  The  mean  and  variance  of  the  intensity  are  given  by 

A  2 

Ug  =  E  [s]  =  2a  ,  (mean) 

3 

(6a) 

3 


(6b) 


m2 


=  E  =  2ij^  ,  (second  moment) 


9  A  9 

Og  =  E  [(s-Us)2]  =  Ug  ,  (variance) 


(6c) 


The  contrast  ratio  (a  measure  of  signal  to  noise)  is  defined  as 


standard  deviation 
mean 


(7) 


For  fully  developed  speckle,  Cj.  =  1  regardless  of  the  signal  strength.  Thus, 
speckle  is  a  signal  dependent  on  multiplicative  noise.  This  is  in  contrast  to 
images  with  additive  noise  where  a  change  in  total  light  level  can  affect  the 
apparent  signal  to  noise  [4]. 

A.  Fully  Developed  Speckle  in  the  Presence  of  Coherent  Background 
(Ricean  Statistics) 

The  effect  of  coherently  adding  a  uniform  intensity  to  the  speckle 
intensity  is  of  more  than  pedantic  interest.  If  the  assumption  that  the  sur¬ 
face  roughness  is  greater  than  a  wavelength  no  longer  can  be  made,  then  the 
resulting  speckle  pattern  can  be  analyzed  by  assuming  that  it  consists  of  a 
uniform  field  coherently  added  to  a  speckle  field  [1,2].  This  would  corre¬ 
spond  to  scattering  from  a  slightly  rough  object  with  the  uniform  field 
corresponding  to  a  specular  reflection  and  the  speckle  field  corresponding  to 
diffuse  scattering.  We  would  expect  to  encounter  this  mixed  beam  case  more 
often  than  the  fully  developed  case  (for  example,  holography). 

In  this  case  the  speckle  complex  amplitude  contains  a  constant  phasor 

so  that 


a  =  as  +  hr  +  j  ai  ,  (8) 

ag  =  constant  with  a  fixed  phase  angle  (zero  phase  angle  can  be  assumed). 

? 

Let  a  =  Sf,.  Then 
s  ^ 


''  7  2 

s  =  (as  +  aR)^  +  a 


=  So  +  s  +  2asaR 

The  statistics  of  s  are  given  as  follows: 


(9) 
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Joint  Density 


p(s, e)  = 


1  I  s  +  Sq  -  2  “y  s  Sq  cos  9  \ 

exp  I - 2 - I.  s>0,  -TT<0<n 

V  / 


0  ,  otherwise. 


Marginal  Density 


s  +  s. 


lo  I  2  \  s  Sq/ Us  I  ,  s  >  0 


0  ,  otherwise. 


Io(.)  =  zero  order  modified  Bessel  function  of  first  kind'. 


Us  ~  Efs]  E[sJ  +  Sq  +  2E 


=  E[sl  +  So  +  2E  [as]  E  [aR] 
-  E[s]  +  Sq  =  Us  Sq  • 


Second  Moment 


E  [(s)2]  =  2u:  +  4usSo  + 


s  ^iJsSo  +  . 


Variance 


2  2 

Os  =  u  (1  +  2r),  where  r  =  beam  ratio  =  Sq/us 


The  contrast  ratio  is  given  by 


Cr  =Jl  +  2r/(l  +  r)  . 


For  large  values  of  r 


B.  N-Look  Averaging 


Cr  =  j2/r 


One  method  used  to  reduce  speckle  noise  is  the  intensity  addition  of 
N  speckle  patterns 


^  SR 

k»l 


5 


If  each  speckle  field  sj^  is  statistically  Independent  of  all  the  other  N-1 
fields,  and  all  the  average  intensities  are  distinct,  i.e.,  4  P£  If 

k  a,  then  the  probability  density  is 

N  N-2 

Uk 

p(s)  =  - exp 

^  (Pk'hp)  «k,p 
p=l 
pj^k 


p(s)  =  0,  otherwise* 

If  the  average  intensity  is  the  same  for  all  N  fields  [1],  i*e.,  pj^  =  Pq  for 
k  =  1 , . • • * ,N,  then 


p(s)  =  - jT  exp 

(N-1)!  PqN 


,  s  >  0 


=  0 


,  s  >  0 


(19) 


The  average  value  of  intensity,  the  variance  and  the  contrast  ratio  may  be 
calculated  from  this  probability  density  function  giving 


<  s  >  *=  Po  (20) 

<  s2  >  =  (1  +  1/N)  po^  (21) 

(22) 

Y  »  1/V^*  (23) 
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III.  DIGITAL  MODEL 


For  our  initial  experiments,  a  very  simple  test  pattern  was  selected. 
This  test  pattern  consists  of  a  64  x  192  array  of  pixels  which  contain  square 
cells  of  sizes  1  x  1 ,  2  x  2,  4  x  4,  8  x  8,  and  16  x  16  at  eight  different  in¬ 
tensities.  This  pattern  is  similar  to  a  chart  developed  by  Rose[8]  f  r  exper¬ 
iments  in  human  vision.  This  type  of  pattern  was  also  used  in  a  previous 
investigation  of  speckle  nolse[9].  Figure  1  shows  the  optical  and  digital 
test  patterns 

If  we  let 

v(x,y)  =  image  intensity  at  x,y 

T(x,y)  =  object  transmittance  (or  reflectivity)  at  x,y 


then 


v(x,y)  =  T(x,y)  s(x,y)  (24) 

is  the  speckle  image  model  in  the  simplest  case  with  no  receiver  noise. 

To  generate  the  speckle  noise  Intensity  s(x,y)  for  fully  developed 
speckle,  we  first  generate  a  pseudo-random  pair  of  Gaussian  random  variables, 
then  calculate  their  magnitude  squared.  This  algorithm,  repeated  for  every 
point  in  the  image  field,  will  produce  a  speckle  noise  field.  One  hundred 
twenty  eight  speckle  noise  fields  64  x  192  pixels  in  size  were  digitally  gene¬ 
rated  . 


In  addition  to  the  noise  fields  with  fully  developed  speckle,  eight 
frames  of  size  64  x  192  that  contain  speckle  noise  with  Rlcean  statistics  were 
generated.  The  beam  ratios  selected  were  0,  2.41,  6.46,  14.48,  30.49,  62.5, 
126.5  and  254.5.  These  beam  ratios  were  selected  to  provide  the  same  beam 
ratio  Cr  as  occurs  in  N-look  averaging  for  N=l,  2,  4,  ....,  128. 

In  the  more  general  case  of  spatially  Invariant  diffraction  limited 
imaging  geometry,  the  observed  image  is  given  by  first  finding  the  amplitude 
field  in  the  imaging  plane  as 


oo  00 

A(x,y)  =  f  /  K(x-j;,  y-n)  a(c,n)d^dn  (25) 

—  00  —00 


where  a(c,n)  is  the  complex  amplitude  distribution  of  the  object.  Then  the 
image  intensity  is  given  by 

v(x,y)  =  |A(x,y)  +  n(x,y)  (26) 

where  n(x,y)  represents  the  additive  detector  noise. 
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IV.  MODEL  VALIDATION:  N-LOOK  FILTER  RESULTS 


Figure  2  shows  the  histogram  of  the  digital  test  pattern  in  the  presence 
of  fully  developed  speckle.  This  histogram  represents  the  convolution  of  the 
test  pattern  Intensity  distribution  with  the  speckle  intensity  distribution 
given  by  (5a).  Figure  3  shows  the  histogram  of  the  Rlcean  speckled  test  pat¬ 
tern  for  Y  =  6.46  (Note  that  y  =  0  would  be  the  same  as  fully  developed 
speckle ) . 

To  validate  the  digital  speckle  models  we  compare  the  N-look,  linear  tem¬ 


poral  averaging  filter  outputs  of  digital  and  optical  experiments.  Let 

vk(x,y)  =  T(x,y)  sk(x,y)  (27) 

where  k  =  1,  ...,  N  and  vj^  represents  the  k^^  image  of  the  object.  Then 

E[vk]  =  TE[sk)  -  Tps  •  (28) 


This  suggests 


TN(x,y) 


JL  1 
Vs  N 


N 

^  vk(x,y) 
k=l 


(29) 


as  an  estimate  of  T(x,y).  This  is  a  point  by  point  (I.e.,  for  each  x,y) 
average  of  N  independent  samples  of  a  random  variable.  Therefore, 


E[Tn]  =  4 


N 

r  E[vk] 

k=l 


T 


(30) 


which  means  Tj^f  is  an  unbiased  estimate  of  T. 


Using  a  probability  model  such 


as  (19)  for  Tf^,  it  is  quite  easy  to  show  that  It  is  the  maximum  likelihood 
estimate  of  T  given  vi^(x,y),  k=l,....,N.  Also, 


E[T 


n1 


1 

2 


x2 

n2 


N 

E  E[SkS£] 
k,)t=l 


(31) 


Since  S]^  and  sjj  are  Independent  from  (6b)  we  get 


(  2 

'  2u^  k  =  £ 


E[skS£]  =  \ 


V 


k  £  . 


(32) 


8 


This  gives 


MT,1  -  Ti  +  t2 


(33) 


From  (33)  we  have 


var[TN] 


E[Tn  -  E(Tn)]2  = 


t2 

N 


(34) 


Therefore  the  contrast  ratio  Is 


Thus  the  contrast  should  improve  by  VlT as  the  number  of  looks,  N  Is  Increased. 

The  above  N  look  digital  filter  was  applied  to  a  sequence  of  128  computer 
simulated  speckle  Images  of  the  digital  test  pattern  of  Figure  1.  Filtered 
outputs  for  N=1 , 2 , 4 , 8 , 16 , 32 , 64  and  128  were  produced  (see  Figure  4).  Output 
images  were  displayed  on  a  high  quality  TV  monitor  and  photographs  were  made 
for  evaluation  by  a  set  of  observers.  The  observers  made  Independent  evalua¬ 
tions  of  the  minimum  detectable  contrast  as  a  function  of  the  object  size  for 
the  digital  Imagery  and  for  optical  Images  produced  (see  Figure  5  [4]  using 
the  test  pattern  In  Figure  1).  Figure  6  [10]  shows  the  comparison  between  the 
optical  and  digital  speckle  experiments  for  minimum  detectable  contrast  for 
objects  whose  size  (area)  is  about  sixteen  times  the  speckle  size.  The 
results  show  agreement  between  the  two  experiments  except  for  the  effects  due 
to  the  optical  transfer  function.  The  reader  can  perform  his  own  evaluation 
by  comparing  the  digitally  produced  data  shown  In  Figure  4  with  the  optically 
produced  data  In  Figure  5.  The  effect  of  the  optical  transfer  function  can  be 
seen  In  the  N-128  case  In  Figure  5  where  the  low  contrast  elements  of  rows  3 
and  4  are  not  visible.  To 'verify  the  Rlcean  speckle  model,  one  needs  to  com¬ 
pare  the  digital  and  optical  Rlcean  speckle  Images  at  various  beam  ratios.  It 
has  been  established  earller[4]  that  a  Rlcean  speckle  with  beam  ratio  r  and  an 
N-look  temporarily  averaged  speckle  have  the  same  contrast  when 

Vl  +  2r  ^ 

1  +  r  Vn 

or  r  =  N  +Vn2  -  N  -  1  .  (36) 

Thus  Rlcean  speckle  at  r  =  0,  2.41,  6-46,  14.48,  30.49,  62.5,  126.5  and  254.5 
corresponds  to  the  N-look  averaged  speckle  for  N  =  1,2,4,8,16,32,64  and  128 
respectively.  Figure  7[10]  shows  the  digitally  and  optically  produced  Rlcean 
speckle  Images  respectively.  An  agreement  consistent  with  Figure  6  was 
obtained  In  this  case  also. 

The  agreement  obtained  between  the  optical  and  digital  Imagery  validates 
the  use  of  the  digital  model  In  speckle  noise  experiments  Involving  human 
observers . 
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Nonlinear  Temporal  Filters 


It  Is  difficult  to  accomplish  nonlinear  filtering  of  speckle  noise  In  an 
optical  system.  However,  In  a  digital  system  it  is  quite  easy  to  produce  non¬ 
linear  filters.  Two  nonlinear  filters  that  were  simple  to  Implement  digitally 
were  selected  for  the  Initial  experlments- 

Square  Root  Filter 


From  Eq.  (27)  we  can  write 


which  suggests 


h  ^V2  V2 

k  k 


tV2  =  e[v^^]/E[s^2] 


TN(x,y)  =  <  7; 


1  ”  V2  V2  I  2 

N  '^k  J  [ 


16  1  ^  V2.  .  ^ 


Squaring  Filter 

Also  from  Eq.  (27) 


2  .,.2  2 
V,  =»  s, 
k  k 


which  gives 


T  =]EtvJ]/E[8^] 


1  I  2 

”  k'l  j 


as  another  estimate  of  T. 
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These  nonlinear  filters  did  not  produce  a  measurable  change  in  the  mini¬ 
mum  detectable  contrast  compared  to  the  linear  temporal  filter.  Observers  did 
comment  however  that  the  noise  in  the  squaring  filter  data  appeared  less 
objectionable.  Experimental  results  using  these  filters  are  shown  together 
with  spatial  averaging  results  of  the  next  section. 


V.  SPATIAL  FILTERS  FOR  SPECKLE  REDUCTION 

The  results  of  the  N-look  temporal  averaging  method  show  that  a  large 
number  of  looks  (N  >  8)  may  be  required  to  detect  objects  which  have  small 
size  and/or  low  contrast.  Spatial  filtering  of  speckle  images  offers  a  means 
of  speckle  noise  reduction  when  the  number  of  available  looks  is  not  very 
large.  We  consider  several  simple  linear  spatial  filters. 

A.  Spatial  Averaging  Filter 

A  standard  technique  used  to  reduce  speckle  noise  in  synthetic  aper¬ 
ture  radar  systems  is  to  Integrate  the  intensity  values  of  adjacent  pixels. 
This  type  of  speckle  filtering  is  digitally  produced  very  easily  in  the 
following  fashion.  Let 


VN(x.y) 


A 

ss 


I 

N 


N 

S  Vk(x,y) 
k=l 


(43) 


and  define  a  spatial  average  Vfj  as 

VN(x,y)  =  V2[vfj(x,y)  +  V4{vn(x  -  A,y) 

+  vn(x  +  A,y)  +  VN(x,y  -  A) 

+  VN(x,y  +  A)}]  .  (44) 

For  digital  images  x,y  are  integer  indices  and  A  =  1.  From  Eqs.  (27)  and  (43) 
we  would  have 


TN(x,y)  =  -p  VN(x,y)  (45) 

Us 

If  T(x,y)  is  almost  constant  over  the  averaging  window.  In  general,  the  size 
of  the  window  will  govern  the  tradeoff  between  resolution  and  speckle  noise 
suppression.  In  these  experiments,  the  spatial  filter  Improved  object  recogni¬ 
tion  for  the  larger  objects  with  an  accompanying  reduction  in  resolution.  The 
spatial  filter  was  combined  with  each  of  the  other  filters  (linear,  square, 
and  square  root)  and  the  results  are  shown  in  Figure  8.  The  minimum  detect¬ 
able  contrast  of  large  objects  for  a  linear  temporal  averaging  of  N  with  spa¬ 
tial  averaging  is  approximately  equal  to  the  same  minimum  detectable  contrast 
for  a  linear  temporal  averaging  of  2N  without  spatial  averaging.  Figure  8 
shows  that,  as  was  the  case  without  spatial  filtering,  the  nonlinear  filters 
produce  no  visible  Improvement. 
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i.  Homomorphic  Filtering  for  Speckle  Reduction 

Since  speckled  images  involve  multiplicative  noise,  linear  filtering 
methods  may  not  be  employed  arbitrarily.  Consider  an  N-look  averaged  speckle 
image  given  as 


VN(x.y)  =  T(x,y)  4(x.y)  ,  T  >  0  ,  S  >  0  (46) 

where  Sjq(x,y)  represents  the  speckle  noise  after  averaging  over  N  looks. 
Taking  the  logarithm  of  both  sides  of  (46)  we  can  write 

WN(x,y)  =  Z(x,y)  +  TiM(x,y)  (47) 

where 


WN(x,y)  =  log  VN(x,y) 
Z(x,y)  =  log  T(x,y) 
nN(x,y)  =  log  SN(x,y) 


Eq.  (47)  now  represents  observations  WN(x,y)  of  a  signal  Z(x,y)  la  the  pre¬ 
sence  of  additive  noise  Ti(x,y).  It  has  been  shown  by  Arsenault  [  11]  that  for 
N  2,  the  additive  noise  may  be  modeled  reasonably  well  by  a  Gaussian  density 
function.  Regardless  of  the  probability  distribution  of  nN>  can  design 

easily  optimum  linear  filters  to  estimate  Z(x,y)  from  Wiij(x,y).  Figure  9  shows 
the  block  diagram  of  the  overall  filter  algorithm.  This  filter  is  also  called 
a  homomorphic  filter.  We  shall  compare  this  filter  with  the  conventional  fil¬ 
tering  method  shown  In  Figure  10.  Here,  we  use  the  same  filter  as  in  homomor¬ 
phic  filtering,  but  logarithmic  and  exponential  operations  are  deleted.  The 
following  filters  are  considered. 


C.  Low  Pass  Filter  (FI) 

An  ad-hoc  digital  low  pass  filter  was  designed  to  filter  the  noise. 
The  frequency  response  of  this  filter  is  given  by 

H(aix,  '*>y)  ~  h(t^)  l'(my)  (^8) 

where  a)x  and  toy  are  the  spatial  frequency  variables  In  the  Fourier  transform 
domain  of  the  sampled  images.  Figure  11  shows  the  normalized  frequency 
response  h(a)).  The  passband  Is  chosen  to  be  [4*4]  and  the  transition 
width  to  be  The  filter  was  Implemented  via  the  discrete  Fourier  trans¬ 

formation  of  the  speckle  Image. 

D.  Wiener  Filter  (F2) 

For  a  signal  with  additive  noise,  the  frequency  response  of  the  opti¬ 
mum  linear  mean  square  filter  (called  the  Wlenei  filter)  Is  given  by 

G(a)x,a)y)  =  S2(a)x»  Wy)/(S2(ux,«y)  +  S^(a)x.Wy))  (49) 
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where  $2  and  Sn  denote  the  spectral  density  functions  of  the  signal  and  noise 
processes.  Since  in  our  digital  simulations  the  speckle  size  was  chosen  to 
coincide  with  the  sampling  grid  size,  the  fully  developed  speckle  field  is 
Independent  from  pixel  to  pixel.  Hence,  nN(*»y)  considered  as  a  white 

noise  field  of  mean  and  \iz,  the  mean  value  of  Z  to  be  zero  by  subtracting 
the  sample  mean  py  of  the  observation  field  from  each  observation  Wfj(x,y). 
Hence,  for  zero  mean  data 


/  7r2 

-  6  •  "  =  ' 


0^,  0^)  “  “  I 


N  •  2 


(50) 


where  we  have  used  the  result  of  Arsenault [ 11]  to  evaluate  a^.  To  estimate 
the  spectral  density  of  the  signal  Z(x,y)  we  assume  its  covariance  function  is 
exponentially  decreasing  in  each  dimension,  l.e., 

Cov[Z(x,y)]  =  E[(Z(x',y')  -  uz)(Z(x+x',y+y')  - 

-  2  hi  hi 

=  Px  Py  (51) 


where  0  £  |  p  |£  1.  For  zero  mean  observation  Wiij(x,y),  the  parameters  py 

can  be  estimated  as 

Px  =  Z  r  W(x,y)  W(x+l,y)/(E  Z(W(x,y))2)  (52a) 

X  y  X  y 

Py  =  E  Z  W(x,y)  W(x,y+1)/(E  E(W(x,y))2)  (52b) 

X  y  X  y 


The  variance  Is  estimated  from  the  sample  variance  of  W(x,y)  as 


2  2  2 

Note  that  this  model  does  not  require  the  a  priori  knowledge  of  the  object 
covariance  since  all  the  measurements  are  made  on  the  available  data.  Experi¬ 
mentally,  this  model  was  found  to  be  satisfactory  although  a  better  covariance 
model  would  be  ultimately  desirable.  The  object  spectral  density  function 
corresponding  to  the  foregoing  covariance  model  is  simply  given  by  its  Fourier 
transform  as 


Sz(<^» 


02  [(1  -  2ax  coscox  ((1  -  2ay  cosujy)]"^ 


(53) 
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where 


r2  =  aj  (1  -  p^)(l  -  pj)/(l  +  pj)(l  +  pj) 

^  X  y  X  y 

2  2 
Ox  =  Px/(1  +  Pjj)  >  ®y  ^  Py/(1  +  Py) 

Given  S2  and  Sy,  the  Wiener  filter  design  Is  complete.  It  Is  implemented  via 
the  discrete  Fourier  transform  of  the  sampled  Image  as  follows: 

1.  Take  the  logarithm  of  the  Image  Vf}(x,y)  to  obtain  Wjj(x,y). 

2.  Find  the  sample  mean  of  WjjCx.y)  and  subtract  from  each  pixel 

value . 

2 

3.  Calculate  p^,  py,  and  the  Wiener  filter  frequency  response. 

4.  Take  the  discrete  Fourier  transform  (DFT)  of  the  Image  (of  appro¬ 
priate  size).  In  the  present  case  (see  examples  In  paragraph  G  below)  a 

256x256  transform  was  adequate. 

5.  Multiply  the  Image  DFT  with  the  Wiener  filter  gain  sampled  at 
appropriate  frequencies. 

6.  Take  Inverse  DFT,  exponentiate  and  scale  the  Image  appropriately 
for  display. 

E.  Spatial  Averaging  Filter  (F3) 

The  homomorphic  spatial  filter  is  implemented  by  spatial  averaging  of 
WN(x,y)  according  to  (44). 

F.  Median  Filter  (F4) 

Here  a  spatial  window  of  suitable  size  (3x3  here)  Is  chosen.  We 
slide  this  window  on  the  Image  field  W{g(x,y)  and  for  each  position,  replace 
the  pixel  value  at  the  center  of  the  window  by  the  median  value  of  all  the 
pixels  lying  In  the  window. 

G.  Examples  and  Comparisons 

The  foregoing  filters  (FI  to  F4)  were  Implemented  digitally  both  In 
the  homomorphic  mode  and  in  the  conventional  mode.  In  the  homomorphic  mode, 
the  logarithmically  transformed  Image  is  the  Input  to  a  spatial  filter  and  its 
output  Is  exponentiated  (Fig  9).  In  the  conventional  mode,  the  logarithmic 
and  the  exponential  transformation  are  deleted.  The  spatial  filters,  once 
designed,  are  kept  the  same  In  either  mode. 

Figure  12  shows  the  test  pattern  used  and  the  corresponding  N-look 
averaged  speckle  images  for  N»l,2,4  that  are  available  for  digital  processing. 
Figures  13  to  16  show  the  results  of  applying  homomorphic  filters  FI  to  F4  to 
these  images.  These  results  show  that  filters  FI,  F2  and  F3  tend  to  smooth 
the  speckle.  The  median  filter  performs  some  smoothing  but  distorts  smaller 
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objects.  In  a  detection  experiment  this  could  cause  high  false  alarm  rates. 
This  Is  not  completely  surprising  because  the  median  filter  Is  known  to  per¬ 
form  well  on  binary  noise  whereas  presently  (In  the  homomorphic  mode)  the 
noise  Is  nearly  Gaussian. 

Figures  17  to  19  compare  the  performances  of  the  various  homomorphic 
filters.  Since  the  noise  Is  additive  and  nearly  Gaussian  as  expected,  we  find 
the  Wiener  filter  to  give  the  best  performance.  The  low  pass  filter  tends  to 
oversmooth  and  somewhat  distort  the  size  and  shape  of  the  objects.  The  median 
filter  affects  most,  the  resolution  of  objects  to  which  are  close  to  the  size 
of  the  median  filter  window.  It  Is  notable  that  both  the  Wiener  and  Spatial 
averaging  filters  are  able  to  resolve  objects  of  smaller  size  (see  the  fourth 
and  the  fifth  rows). 

Finally,  we  compare  the  conventional  vs  homomorphic  filters  In 
Figures  20  to  23.  For  larger  objects  (low  resolution),  the  performance  of  the 
linear  and  homomorphic  filters  Is  nearly  equivalent.  However  for  smaller  size 
objects  (rows  3,  4,  and  5)  l.e.,  at  higher  resolution,  the  homomorphic  filters 
definitely  seem  to  perform  better.  This  Is  particularly  noticeable  for  the 
Wiener  and  Spatial  averaging  filters  (Figures  21,  22). 

H.  Conclusion 

Based  on  the  foregoing  experiments  we  conclude  the  following: 

1.  Homomorphic  filtering  has  a  definite  advantage  over  straightfor¬ 
ward  linear  filtering  for  speckle  reduction. 

2.  The  homomorphic  Wiener  filter  yields  most  favorable  results  In 
terms  of  resolution  vs  noise  reduction.  More  accurate  object  spectral  density 
models  could  be  useful  In  the  design  of  Wiener  filters. 
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Figure  1.  Optical  and  digital  test  patterns. 
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(b)  DIGITAL  TEST  PATTERN 
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Figure  2,  Histogram  of  the  digital  test  pattern 
with  fully  developed  speckle  noise. 


Figure  4. 
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Figure  6.  Contrast  required  to  detect  objects 

sixteen  times  the  speckle  size.  [10] 
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Figure  7.  Optical  and  digital  simulations  of  Ricean  speckle  images.  (Page 


Figure  7.  Optical  and  digital  simulations  of  Ricean  speckle  Images.  (Page  2  of  2) 
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Figure  16.  Homomorphic  median  filter, 


FILTERING  AFTER  LOGARITHMIC  TRANSFORMATION 
FULLY-DEVELOPED  SPECKLE  IMAGE 


Figure  17.  Comparison  between  various  homomorphic  filters. 


FILTERING  AFTER  LOGARITHMIC  TRANSFORMATION 

2-LOOK  AVERAGED  IMAGE 
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Figure  18.  Comparison  between  various  homomorphic  filters. 


FILTERtIMG  AFTER  LOGARITHMIC  TRANSFORMATION 
4-LOOK  AVERAGED  IMAGE 
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SPATIAL  AVERAGE  FILTER 


> 


Figure  21.  Comparison  of  linear  and  homomorphic  spatial  averaging 
filters,  N  =  1 ,2, A. 


WIENER  FILTER 


INEAR  INPUT  LOG  INPUT 

‘6  22.  Comparison  of  linear  and  homomorphic  Wiener  filt 


MEDIAN  FILTER 


LINEAR  INPUT  LOG  INPUT 

Figure  23.  Comparison  of  linear  and  homomorphic  median  filters. 
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